Unified dark energy models : a phenomenological approach 
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A phenomenological approach is proposed to the problem of universe accelerated expansion and 
of the dark energy nature. A general class of models is introduced whose energy density depends 
on the redshift z in such a way that a smooth transition among the three main phases of the 
universe evolution (radiation era, matter domination, asymptotical de Sitter state) is naturally 
achieved. We use the estimated age of the universe, the Hubble diagram of Type la Supernovae and 
the angular size - redshift relation for compact and ultracompact radio structures to test whether 
the model is in agreement with astrophysical observation and to constrain its main parameters. 
Although phenomenologically motivated, the model may be straightforwardly interpreted as a two 
fluids scenario in which the quintessence is generated by a suitably chosen scalar field potential. On 
the other hand, the same model may also be read in the context of unified dark energy models or 
in the framework of modified Friedmann equation theories. 
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I. INTRODUCTION 

In the last few years, an increasing bulk of data has 
been accumulated leading to the emergence of a new 
cosmological scenario. The Hubble diagram of type la 
Supernovae (SNela) first indicated that the universe ex- 
pansion is today accelerating Q, Q- The precise de- 
termination of first and second peaks in the anisotropy 
spectrum of cosmic microwave background radiation 
(CMBR) by the BOOMERanG and MAXIMA collab- 
orations 0] strongly suggested that the geometry of the 
universe is spatially flat. When combined with the data 
on the matter density parameter Qm, these results lead 
to the conclusion that the contribution Qx of dark en- 
ergy is the dominant one, being Q,m — 0.3, Ox — 0.7. 
This picture has been strenghtened by the recent deter- 
mination of CMBR spectrum measured by the WMAP 
team 

According to the standard recipe, pressureless cold 
dark matter and a homogenously distributed cosmic fluid 
with negative pressure, referred to as dark energy, fill the 
universe making up of order 95% of its energy budget. 
What is the nature of this dark energy still remains an 
open and fascinating problem. The simplest explanation 
claims for the cosmological constant A thus leading to 
the so called ACDM model 5]. Although being the best 
fit to most of the available astrophysical data 0, the 
ACDM model is also plagued by many problems on dif- 
ferent scales. If interpreted as vacuum energy, A is up 
to 120 orders of magnitudes smaller than the predicted 
value. Furthermore, one should also solve the coincidece 
problem, i.e. the nearly equivalence of the matter and A 
contribution to the total energy density. 

As a response to these problems, much interest has 
been devoted to models with dynamical vacuum energy, 
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dubbed quintessence 6]. These models typically involve 
scalar fields with a particular class of potentials, allow- 
ing the vacuum energy to become dominant only re- 
cently (see 0, for comprehensive reviews). Altough 
quintessence by a scalar field is the most studied can- 
didate for dark energy, it generally does not avoid ad 
hoc fine tuning to solve the coincidence problem. On 
the other hand, a quintessential behaviour may also be 
recovered without the need of scalar fields, but simply 
by taking into account the effective contribution to cos- 
mology of some (usually neglected aspects) of fundamen- 
tal physics ;9j. A first tentative were undertaken show- 
ing that a universe with a non vanishing torsion field 
is consistent with SNela Hubble diagram and Sunyaev- 
Zel'dovich data on clusters of galaxies ^(j- The same 
quintessential framework can be obtained with the ex- 
tension of Einstein gravity to higher order curvature in- 
variants leading to a model which is in good agreement 
with the SNela Hubble diagram and the estimated age 
of the universe • ^ is worth noting that these alter- 
native schemes provide naturally a cosmological compo- 
nent with negative pressure whose origin is simply related 
to the geometry of the universe itself thus overcoming 
the problems linked to the physical significance of scalar 
fields. 

Despite the broad interest in dark matter and dark en- 
ergy, their physical properties are still poorly understood 
at a fundamental level and, indeed, it has never been 
shown that the two are in fact two different ingredients. 
This observation motivated the great interest recently de- 
voted to a completely different approach to quintessence. 
Rather than the fine tuning of a scalar field potential, it 
is also possible to explain the acceleration of the universe 
by introducing a cosmic fluid with an exotic equation of 
state causing it to act like dark matter at high density 
and dark energy at low density. An attractive feature of 
these models is that they can explain both dark energy 
and dark matter with a single component (thus automati- 
cally solving the coincidence problem) and have therefore 
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been referred to as unified dark energy (UDE) or unified 
dark matter (UDM). Some interesting examples are the 
generalized Chaplygin gas [l^, the tachyonic field [lflj 
and the condensate cosmology [l4f . 

It is worth noting that all the dark energy models (both 
with scalar fields or UDE) proposed up to now predict 
that the expansion of the universe is a two phase pro- 
cess : it is first determined by a matter - like term and 
it is then driven by the quintessence - like component to- 
wards an asymptotically de Sitter state. However, it is 
well established that there is also a third phase preced- 
ing these two, i.e. the radiation dominated era. It is thus 
interesting to look for a model which is able to predict a 
smooth transition from one phase to the following one in 
a natural way. 

A quite confused picture emerges from the previous 
discussion about dark energy and its nature. Many mod- 
els have been proposed to succesfully reproduce the astro- 
physical observations available up to date, but they are so 
different each other that the mistery of the dark energy is 
far to be solved. Given the state of the art, a different ap- 
proach to the problem is welcome. We think that a first 
step toward understanding the nature of dark side (dark 
energy and dark matter) of the universe is to explore phe- 
nomenological models which are able to reproduce what 
we observe. To this aim, we consider a single fluid whose 
energy density scales with the redshift in such a way that 
the radiation dominated era, the matter domination and 
the accelerating phase have been naturally achieved. We 
constrain the model parameters using the estimated age 
of the universe, the SNela Hubble diagram and the an- 
gular size - redshift relation of radio structures. Although 
phenomenologically motivated, nonetheless the model we 
propose may be physically interpreted in terms of an ef- 
fective two fluids scenario with the dark energy compo- 
nent represented by a scalar field with a suitably chosen 
interaction potential. On the other hand, the models may 
also be considered as a new member of the UDE class. 

The plan of the paper is as follows. In Sect. II we 
present the main feature of the class of phenomenological 
models we propose. Some considerations on the allowed 
range for the model parameters lead us to restrict our 
attention (at least in this first analysis) to a particular 
class of models for the reasons we explain in Sect. III. The 
estimated age of the universe, the SNela Hubble diagram 
and the angular size - redshift relation for compact radio 
structures are used in Sects. IV, V and VI to test whether 
the model is a viable one and to constrain its parameters. 
Sect. VII is devoted to the interpretation of the model in 
terms of an effective two fluids model with a scalar field 
whose interaction potential is reconstructed. Finally, we 
summarize and conclude in Sect. VIII. 



II. A GENERAL CLASS OF MODELS 

A phenomenological approach is proposed here to build 
a model which is able to fit the available data and leads 



to an accelerated expansion. To this aim, let us observe 
that most of the cosmological models (both with one or 
two fluids) predict that, during its evolution, the universe 
undergoes first a radiation dominated expansion (i.e. the 
energy density p scales with the scale factor R as i?~ 4 ), 
then a matter dominated phase (with p ~ _R~ 3 ) and fi- 
nally a de Sitter - like expansion with the energy density 
asymptotically approaching a constant value. This con- 
sideration leads us to assume the following expression for 
the energy density 1 which our model is made of : 



p(R) = A{l + ^ 











1 + 











(1) 



with < a < f3, s and b (with s < b) two scaling factors 
and A a normalization constant. For several applications, 
it is useful to rewrite the energy density as a function of 
the redshift z. Replacing R = (1 + z)^ 1 in Eq.QJ, we 
get: 



p(z)=A 1 + 



having defined : 



1 + z 

1 + Z S 



0-a r 



1 + z 

l + Zb 



(2) 



(3) 
(4) 



z s = 1/s - 1 , 

z b = 1/b - 1 . 
It is quite easy to see that : 

p ~ R-P for R«s , 
p ~ R~ a for s « R « b , 



p ~ const for R >> b . 

Choosing (a, j3) = (3,4), the model we obtain is able to 
mimic a universe undergoing first a radiation dominated 
era (for z » z s ), then a matter dominated phase (for 
Zb « z « z s ) and finally approaching a de Sitter phase 
with constant energy. This is just what we need. How- 
ever, for sake of completeness, we discuss in this section 



1 As yet said in the introduction, the model we are considering 
may be considered as composed by one or two fluids. Here, 
we prefer the single fluid interpretation, but the results we will 
obtain hold whatever is the number of fluids the model is made of. 
In Sect. VII, we will investigate the consequences of our results on 
the scalar field potential, if the model is considered as composed 
by two fluids. 
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the main properties of the model for the general case, i.e. 
with (a,/3) not fixed 2 . 

As a preliminary step, let us recall the Friedmann equa- 
tions 



8ttG 



R 9 k 
2 R +H+ R? 



SivGp 



(5) 



(6) 



where H = R/R is the Hubble parameter, p the pressure 
and the dot denoting the derivative with respect to t. 
From now on, we will consider only flat models so that 
k = 0. Inserting Eq.QJ into Eq.© for a flat universe 
and evaluating it today (where we set R = 1), we may 
express the normalization constant A as : 



A 



Pcrit 



{l+s)P- a {l + b a ) 



(J) 



with Hq the present day Hubble constant and p cr it = 
3Hq /8ttG the critical density. From now on, today eval- 
uated quantities will be denoted by the label "o". Also 
the continuity equation : 



p + 3H(p + p) = 



(8) 



has to be taken into account. 

Using the obvious relation dp/dt = dp/dR x dR/dt and 
the defmtion of the Hubble parameter, it is immediate to 
get the following expression for the pressure : 
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qo 



FIG. 1: Zb as function of qo for the model with (a, f3) = (3, 4) 
and z s = 3454. 



with w = pj p as usual. The expression, as a function of 
z, may be easily obtained replacing R = (1 + z)^ 1 . The 
result for the case (a, f3) = (3, 4) is : 



w(z, a = 3,/3 = 4) 



\l+z b J 



1+.-. 



3 1 



(ii) 

Note that the barotropic factor w strongly depends on 
the redshift z. In particular, Ea. (|ll|l shows that : 

w ~ 1/3 for z >> z s , 



(9) 



which holds whatever is the cosmological model 3 pro- 
vided that H is not vanishing everywhere (i.e. the uni- 
verse is not stationary). Inserting Eq.JTJ into Eq.JjJJ, 
after some algebra we get : 

_ [(a - 3)R +{13- i)s]b a - [3(R + s) + (a - f3)s]R a 



3(R + s)(R a + b a ) 



(10) 



We may refer to this class of models as the Hobbit models. In 
Tolkien's trilogy "The Lord of the Rings", the Hobbits look like 
a "mixture" of the three main people of the book having the 
aspect of Men, (almost) the height of Dwarfs and pointed ears 
as Elfs. In the same way, our models behave as the three main 
fluids of the standard cosmological model. 

Actually, if the model contains more than a single fluid, this 
relation strictly holds only for the total pressure and the total 
energy density. However, it is still valid for each fluid if each 
one satisfies the continuity equation separately. This happens 
whenever the fluids are not interacting as in many quintessence 
models. 



W ~ for Zft << Z << Z s 



w ~ — 1 for z << Zb • 



Before discussing in more detail the behaviour of w with 
z, let us first derive the expression of the deceleration 
parameter q. Combining the two Friedmann equations 
for a fiat case, we easily get : 



da 1 3 p 
a z 2 2 p 

Inserting Eq. i|10|l into this relation gives : 



(12) 



[(a - 2)R +((3- 2)s]b a - [2{R + s) + (a - (3)s}R a 



H 2(R + s)(R a +b a ) 

inserting R = 1 gives the present day value as : 



(13) 



90 



(y - l)a + z s [ay~ 2(1 + y)] + (fl - 4)(1 + y) 



2(2 + z s )(l + y) 



(14) 
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FIG. 2: The barotropic factor w as function of log(l + z) for 
the model with (a,/3) = (3,4) and (qo,z s ) = (-0.5,3454). 

with y = (l+Zb)~ a . It is convenient to solve Ea. (|14f> with 
respect to Zb in order to express this one as a function of 
qo and z s . It is : 



Zb 



a(l + z s ) + (3 - (2 + z s )(2q Q + 2) 



l/c 



1 . (15) 



a - (3 + (2 + z s )(2q + 2) 

In Fig.m wc plot Zfc(qo) for the model with (a, (3) — (3, 4) 
having fixed z s — 3454 (see later for the motivation of 
this choice). For q ~ —0.5, it is z& ~ 0.26 so that such 
a model describes a universe which is dominated by a 
radiation- like fluid for z » 3454, then its dynamical 
evolution is driven by a matter -like fluid until z ~ 0.26 
when a term similar to the cosmological constant begins 
to dominate leading asymptotically to a de Sitter phase. 
This behaviour is quite similar to what is predicted by 
a cosmological model with dark matter and quintessence 
thus showing that the model we are considering is phe- 
nomenologically equivalent to the standard framework. 
Similar considerations refer to the models with other val- 
ues of the parameters (a, (3). 

Some straightforward physical considerations allow us 
to use Ea.l|T5)l to derive constraints on q . First, we note 
that Zb is obviously a finite quantity so that we have to 
reject all values of go which makes the right hand side of 
Ea. (|15|l to diverge. Secondly, it is reasonable to assume 
that Zb > since we need a matter dominated universe in 
the past so that structures can efficiently form. Impos- 
ing these two constraints, Eq. l|15|) leads to the following 
condition : 



QO.min < 50 < 10,- 



with : 



go,-, 



(3-a 



go,-, 



2 + z s 

az s + 2(3 
2(2 + z s ) 



(16) 

(17) 
(18) 



FIG. 3: Derivative of the barotropic factor w(z) vs log (1 + z) 
for the model with (a, fS) = (3, 4) and (q , z s ) = (-0.5, 3454). 
We consider four redshift ranges : (0, 5) (top left), (5, 15) (top 
right), (15,1500) (bottom left), (1500,4500) (bottom right). 



It is possible to see that both go, mm and qo, ma x are almost 
independent on z s for z s G (1000,4000) for the model 
with (a, (3) = (3,4) so that, with a very good approxi- 
mation (more than 0.1%), we may fix: 

qo,mm(a = 3,/3 = 4) ~ -1 , 

<Zo, m az(a = 3,/3 = 4)~-0.25. 

Before constraining the model parameters with the 
available observations, we turn back again to the 
barotropic factor w in order to qualitatively discuss some 
its interesting features. To this aim, we fix (a, (3) = (3, 4) 
and (q a ,z s ) = (-0.5, 3454) giving z b = 0.26 and, in Fig.[H 
we plot w vs log(l + z). Note that w(z) starts from a 
value near to that of a radiation dominated universe (i.e., 
w = 1/3), but it is exactly equal to that value only for 
z ~ lOOzs thus suggesting that the radiation dominated 
phase of the evolution takes place only at the very be- 
ginning. For z < z s , w is almost null (being w = for 
z ^ 11) coherently with the picture of a matter domi- 
nated universe. For z < Zb, w quickly declines towards 
the asymptotic value w = — 1 so that the pressure be- 
comes negative and the universe enters a phase of accel- 
erated expansion. For the chosen values of parameters 
(a, (3) and (qo,z s ), it is w(z = 0) = —2/3. It is also in- 
teresting to look at the derivative of w{z) that we plot in 
Fig-El The derivative is almost vanishing for most of the 
past history of the universe, being significantly different 
from zero only in the recent past because of the transition 
of the cosmic fluid from matter - like to quintessence - like 
regimes. This should suggest the use of an approximated 
equation of state as p — w e ffp with w e ff evaluated as 
a mean for w(z) over the redshift. However, this proce- 
dure must be avoided since it could lead to serious sys- 
tematic errors. Actually, looking at the relative variation 
w^ 1 x dw/dz shows that w strongly depends on z so that 
the introduction of a constant effective w has no physical 
justification. 
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A. Modified Friedmann equations ? 

Up to now, we have derived the main properties of the 
model by implicitely assuming that Eq.JIJ describes the 
energy density of a single fluid accounting for both dark 
matter and dark energy. This phenomenologically moti- 
vated assumption may also be abandoned in favour of a 
different approach to the cosmic acceleration. As recently 
suggested, one might also think that, in the words of [lfij. 
the observed acceleration is not the manifestation of yet 
another new ingredient in the cosmic gas tank, but rather 
a signal of our first real lack of understanding of gravi- 
tational physics. In this framework, one assumes that 
standard matter is the only component of a flat universe, 
while the Friedmann equation |JSJ is replaced by : 



x + (1 — il m .o)- Also limiting to the case (a, (3) = (3, 4), 
our model represents a generalization of the cosmologi- 
cal constant scenario since Ea. (l22|) reduces to this very 
special case only for s = and adjusting b. 

In the following, we still prefer to interpret our model 
as a phenomenological one in the framework of a unified 
description of dark matter and dark energy. However, it 
is worth stressing that, since the astrophysical tests we 
will consider later are mainly sensitive to the shape of the 
Hubble function H(z), the main results we will obtain are 
independent on what is the preferred physical meaning of 
the model among the different possibilities (UDE, matter 
plus scalar field or modified Friedmann equations). 



H 2 = H 2 g(x) 



(19) 



with x = Pm/pcrit and p m scaling as usual (i.e., p m °c 
i?~ 3 ). The function g(x) reduces to x in the early stage 
of the universe evolution, while takes a different (nonlin- 
ear) form later. By suitably chosen g(x), different mod- 
els fitting the astrophysical data may be obtained, the 
most interesting ones being the ( gene ralized) Cardassian 
model [l7| and the DGP gravity (Isj. 

Our phenomenological model may be interpreted in 
this framework provided that we choose : 



(x/n mfi )p/ 3 

(1 + s)P- a (l + b a ) 



n 



m,0 



-1/3- 



"m,0 



-a/3 



(20) 



having used Eq.Q to fix the normalization constant A 
and being fl m , = p m ,o/Pcrit- In particular, for (a,/?) = 
(3,4), Eq.® reduces to: 



^m,0(l+s)(l + 6 3 ) 



+ 9nl(x) (21) 



with 



9m{x) = 



l + s(a;/n m ,o) 1 / 3 (l + 6 3 a;/n m ,o) 
(l + s)(l + 6 3 ) 



(22) 



Note that g(x) ~ x in the early universe as expected 4 . 
It is also interesting to compare Eqs.(j2U and (1221) with 
the corresponding expression for a two fluids model com- 
posed of matter and cosmological constant having g{x) oc 



4 This is easy to show in the case (a, (3) = (3, 4) remembering that 
(s/fi m ,o) 1/3 = 1/R = 1 + Z. 



III. THE SPACE OF PARAMETERS 



The general expression of energy density we are consid- 
ering is characterized by five parameters which we may 
choose to be the two slopes (a,/3), the scaling redshift 
z s , the present day value of the deceleration parameter 
<7o and the Hubble constant Hq entering through the nor- 
malization coefficient A in Eq. . The astrophysical ob- 
servations available up to date could be used to constrain 
the model parameters, but this is a daunting task given 
the large space of parameters to be search for. Actu- 
ally, we have seen that models with (a,/3) = (3,4) are 
the most interesting ones since the energy density scales 
with R as in the standard cosmological framework with 
matter and dark energy. Hence, hereinafter we will fix 
(a,/?) = (3,4). The space of parameters to search for 
is now significantly reduced since we have to consider 
only three of five quantities, namely (z s ,qo, H ). Actu- 
ally, z s marks the transiton of the fluid from radiation - 
like to matter -like regimes. In the standard framework, 
radiation and matter give the same contribution to the 
energy budget at z eq so that it is reasonable to assume 
that our fluid behaves as radiation much before this era. 
Henceforth, a possible choice could be z s ~ z eq . Fitting 
the ACDM model to the CMBR anisotropy spectrum, 
the WMAP collaboration has found (as best fit values) 
z eq = 3454 [4| so that we fix z s to this value. However, 
we have checked that varying z s in the range (1000, 5000) 
does not change the constraints on (q ,H ) being all the 
tests we will discuss later fully degenerate with respect 
to z s . This is not an unexpected result since the avail- 
able observations probe a redshift range which is very far 
from z s whatever its exact value is. These considerations 
leave us with only two unknown quantities to constrain : 
the deceleration parameter qg and the Hubble constant 
Ho. In the following sections, we will use different astro- 
physical observations to check whether the model may be 
reconciled with them and to constrain these two param- 
eters. 
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IV. THE AGE OF THE UNIVERSE 

In order to narrow the parameter space (qo,Ho), as 
a first step, we may compare the predicted age of the 
universe with the constraints coming from both astro- 
physical estimates and WMAP data. Inserting Eq.QJ 
into Eq.|J5J| with A given by Eq.Q, we get : 



to 



l (l + s)' 3 - a (l + b a ) 



1 R- f} {R + s y- a {R a + b a ) 



R 



dR . (23) 



Ea. (|23[l is not analitically solvable, but may be easily 
integrated numerically provided that the values of pa- 
rameters (a,0,z s ,qo,Ho) have been given. Having yet 
fixed (a, (3, z s ) = (3,4,3454), we have to choose only 
the ranges for the two remaining unknown quantities, 
(go j Ho). We let qo vary in the full range determined be- 
fore imposing the physically motivated constraint z& > 0, 
i.e. we take go G ( — 1,-0.25), while we examine mod- 
els with Hq G (56,88) km s _1 Mpc -1 since this is the 
2cr confidence range determined by the final result of 
the HST Key Project 0. Ea.{25)l is integrated over 
a grid in the (qo,Ho) plane with steps of 0.01 in go and 
0.5 in Hq and interpolated for other values. In Fig.01 
we plot the age contours in the (Ho,qo) plane. Super- 
imposed, we show also the contours relative to the 3tr 
confidence range determined by the WMAP data giving 
to G (13.1,14.3) Gyr, in good agreement with other in- 
dependent astrophysical estimates. 

Fig.0] shows that the age test is unable to put any 
constraint on the Hubble constant Hq, while put only a 
lower limit on g imposing the constraint go > —0.85. On 
the other hand, the test is a first evidence that the model 
we are considering is a reliable one since it predicts an 
age of the universe in agreement with the observational 
constraints. 



V. THE SNEIA HUBBLE DIAGRAM 

In order to further constrain the parameters (qo,Ho), 
we fit the model to the Hubble diagram of Type la Su- 
pernovae using the data recently released by the High - z 
Team [2^ and the If A Deep Survey 21] . Following the 
method described in |2dj] , we minimize the quantity : 



1 



N -2 

N 



i=l 



(log H d)i - \ogcd L (zi) + log/i, 



li5 



(24) 



with (\ogHod)i the measured value of the distance (av- 
eraged over the different method used to compute it by 





\\ \\ \ \ 




-0.3 


i\\ V \ \ 




-0.4 


\ A. \ \ \ 

\ \ *\ \ \ 






\ ^ \ \ \ 




-0.5 


\ ^\ \ \ 




o 


\\ v \ 




tr 

-0.6 


\\\ \ \ 










-0.7 








>v VS. 










-0.8 













FIG. 4: Age contours in the (Ho, go) plane for the model with 
(a, (3) = (3,4) and z s = 3454. The age values run from 11 
(the upper line) to 15 Gyr (the lowest one) in steps of 1 Gyr. 
The dashed lines mark the region with to G (13.1, 14.3) Gyr. 



the High-z Team), dh the dimensionless luminosity dis- 
tance, c the speed of light, the Hubble constant in 
units of 65 km s^ 1 Mpc" 1 and <n the reported error and 
the sum is over the TV SNela observed. For our model 
with (a, (3, z s ) = (3, 4, 3454), it is : 



d L (z) = (1 + z) y/(l + s) (1 + b 3 ) x 



K Jo dC 1 1 



i + C 



-1/2 



i + C 

1 + Z b 



-1/2 

(25) 



Note that d^ depends on go through Zf, since this is eval- 
uated using Eci. (|15J) . We fit the model to the data using 
a sample of 162 SNela, 130 from Tonry et al. H3 and 23 
from Barris et al. |2l]]- These have been selected from a 
larger sample according to the two criteria z > 0.01 and 
Ay < 0.5 (being Ay the absorption in the V band) as in 
|20j . Note that we do not use the 42 SNela observed by 
the Supernova Cosmology Project (SCP) since their 
distance modulus have been estimated using a completely 
different approach with respect to those implemented by 
the High - z Team and also used by the If A Deep Survey. 
Actually, Tonry et al. [2(J have shown that inclusion of 
the SCP SNela does not alter the main results of fit- 
ting a model to the Hubble diagram so that we prefer 
to work with a homeogenous dataset even if this lowers 
the sample. We also take care of velocity uncertainties 
in estimating the supernova redshift adding 500 km/s di- 
vided by the redshift in quadrature to following the 
prescription in 20] . The main results of the fitting pro- 
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FIG. 5: The 68 and 95% confidence ranges in the (Ho,qo) 
plane from fitting the model with (a, ft, z s ) = (3,4,3454) to 
the SNela Hubble diagram. 

cedure are presented in Fig.[5] where we plot the 68 and 
95% confidence levels in the (H ,qo) plane. 

Fig.0 shows that fitting to the SNela Hubble diagram 
does not allow us to put constraints on q since the con- 
tours are not closed along the qo direction. This is not 
true for H so that we may get an estimate of this param- 
eter. To this aim, since x 2 is °f order 1, we first define 
the likelihood function as : 



L H HH L L max LqHqoL L max 




FIG. 6: Marginalized likelihood functions (normalized to their 
maximum values) from fitting the model with (a,/3, z s ) = 
(3, 4, 3454) to the SNela Hubble diagram. 

H e (53.1,80.1) km s" 1 Mpc" 1 (95%CL) . 

This is in good agreement with the values obtained by 
the HST Key Project using various standard candles 
and from other independent techniques such as time de- 
lays in multiply imaged lens systems [22] and Sunyaev- 
Zeldovich clusters 231. 



VI. THE ANGULAR SIZE- REDSHIFT TEST 

The relation between the (apparent) angular size and 
the redshift for compact radio structures in quasars and 
radio galaxies has been recently proposed as a possible 
cosmological test. Since radio data probe the redshift 
range bewteen z = 0.011 and z = 4.72, it is clear that 
this test is potentially able to discriminate among differ- 
ent cosmological models breaking the degeneracy present 
in the lower redshift range proven by the SNela Hubble 
diagram. To see how this test works, let us consider an 
object having an intrinsic linear size I and let z be its 
redshift. The apparent angular size is (see, e.g., [24j) : 



C(q ,H ) cxexp{- X 2 /2} (26) 
and then marginalize to get the two following functions : 



£ q (qo) °c J C(q a ,H a ) dH , 
Ch(Ho) oc J C(q a ,H a ) dq ■ 



(27) 



(28) 



In Fig.El we plot the marginalized likelihoods normalized 
to their maximum values. From this plot, we see that the 
SNela test is completely degenerate with respect to qo. 
The maximum is attained for qo = —0.42, but C q varies 
less than 10% over the full range (—1,-0.25). On the 
other hand, Ch strongly depends on Hq (with the maxi- 
mum obtained for Hq = 64.3 km s _1 Mpc -1 ) so that we 
get the following estimates for the Hubble constant : 



8(z) 



I 



l{l + zf 



D A {z) D L {z) 
Ic {l + zf D(l + z) 2 



Hq d L (z) 



dh{z) 



(29) 



with Da(z) the angular diameter distance, D — Ic/Hq, 
the intrinsic angular size (in mas) and we have used the 
relation Dl = (1 + z) 2 Da, with Dl = c/Hq the lu- 
minosity distance yet introduced in Sect. 5. It is worth 
stressing that Eq. H29JI implicitely assumes that the intrin- 
sic linear size I may be considered as a standard rod, i.e. 
it is the same whatever are the properties of the radio 
source. Actually, the validity of this hypothesis is still 
to be demonstrated, both observationally and on a theo- 
eretical ground, and one should consider also a possible 
dependence of I on the the total luminosity L and/or on 
the redshift z. A simple way to parametrize these effects 
is to replace Eq. (2.5]) with the following phenomenological 
one [23 : 



H Q € (58.8,72.3) km s" 1 Mpc" 1 (68%CL) , 



9(z) = K D V< (l + z) n D^\z) (30) 
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FIG. 7: The 68, 95 and 99% confidence ranges in the (D, q Q ) 
plane from fitting the model with (a,/3,z s ) = (3,4,3454) to 
the angular size - redshift relation for compact radio sources. 



with k a normalization constant and (7, n) unknown pa- 
rameters. Gurvits et al. [24j have fitted this relation 
(assuming a A = cosmological model) to the data com- 
ing from 145 radio sources (smoothed in 18 redshift bins) 
selected according to some selection criteria (see later) 
from a sample of 330 5 GHz VLBI sources. Their analysis 
shows that the estimates of the cosmological parameters 
obtained for different choices of (7, n) are consistent with 
each other within the errors. It is worth noting, however, 
that the uncertainties are quite large so that their result 
should be considered as an evidence of the degeneracy 
among different values of (7, n), not as a probe of the an- 
gular size - redshift test being independent on the choice 
of (7,^1). Nonetheless, many authors 25] usually assume 
(7, n) = (0,0) so that we follow this approach and con- 
sider as unknown only the intrinsic angular size D and 
the model parameters. In particular, being this test in- 
dependent on the Hubble constant Hq (since it has been 
included in the D quantity) and having fixed as before 
(a, (3, z s ), the deceleration parameter go is the only model 
parameter we may constrain. 

Following [24[, we only select from the sample in 
Gurvits et al. the radio sources that have spectral in- 
dex a s £ (—0.38,0.18) in order to reduce the intrinsic 
scatter in the angular size - redshift relation and smooth 
the data in (nearly) equally populated redshift bins. We 
minimize the quantity : 
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FIG. 8: Marginalized likelihood functions (normalized to their 
maximum values) from fitting the model with {a, beta, z a ) = 
(3,4,3454) to the angular size - redshift relation of compact 
radio sources. 



with 8 bs,i the observed value of the angular size in the 
redshift bin z i: 9 mo d{zi) given by Eas-PS]! and i(29"|) and 
<7i the uncertainty. Since the resulting \ 2 for the best 
fit models is much lower than 1, we renormalize the er- 
rors in such a way that, for the best fit model, it is 
X 2 = 1. Although not statistically correct, this is the 
usual approach followed when dealing with data affected 
by likely overestimated errors. Moreover, this allows us 
to define marginalized likelihood functions proceeding in 
the same way as with Eas. (|27|) and i|28|) in Sect. 5. The 
main results of the angular size - redshift test are resumed 
in Figs.H and |H1 showing, respectively, the 68%, 95% and 
99% confidence ranges in the (D, g ) plane and the two 
marginalized likelihood functions (normalized to their 
maximum values). Because of the large errors affecting 
the data, we are still not able to constrain the value of 
the deceleration parameter. The best fit value turns out 
to be go = —0.88, but the shape of the likelihood function 
only allows to put an upper limit, go < —0.52 at the 68% 
confidence limit. On the other hand, the best fit value 
of the intrinsic angular size is D — 0.81 mas and we get 
the following constraints : 



(31) 



D e (0.44, 1.48) mas , 



D e (0.31,2.00) mas . 

To investigate the effect of possible systematic errors and 
selection effects, we repeat the angular size - redshift test 
using a different sample comprising only ultracompact 
radio sources given by Jackson |26|. According to him, 
this dataset is more homogeneous than the one in and 
has also been corrected for any selection effect. Jackson 
also gives the values of (z, 9) to be used in the angular 
size - redshift test, while, following his prescription, the 
error on each of the six data points is estimated so that 
X 2 = 1 for the best fit model. We give the correspond- 
ing marginalized likelihood functions defined before in 
Fig-El It is remarkable that the two likelihood functions 
are quite narrow so that it is possible to get constraints 
on both D and go- The best fit value for the deceleration 
parameter is go = —0.64, while the 68% and 95% turn 
out to be : 

go S (-0.76,-0.54) , 
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FIG. 9: Same as Fig.|S]but using the data for ultracompact 
radio sources in 261. 



qo e (-0.83,-0.37) . 

These ranges do not contradict the results obtained be- 
fore from the same test with different data. They are, 
however, significantly narrower. This is in line with what 
is found in ref. [26| for the ACDM model. It seems that 
the sample provided by Jackson allows to narrow the con- 
traints on the fitting parameters (both for our model and 
the ACDM model) because of the lower dispersion of the 
data due to the removal of selection effects. 

Regarding the intrinsic angular size, we find D = 
0.28 mas as best fit value, while the corresponding 68% 
and 95% ranges are: 



px/pi = (b/RY 



for R»b 



D e (0.22,0.34) 



D e (0.15,0.39) mas . 

The intrinsic angular size turns out to be much smaller 
than what has been obtained using the sample in |24| . 
However, this could be the result of the different selection 
criteria used to build the sample (ultracompact instead 
of compact sources). Furthermore, the apparent angular 
size is defined differently for the two samples (see the 
final remark in |2||). 



VII. RECONSTRUCTION OF THE SCALAR 
FIELD POTENTIAL 

The model we have described and tested against some 
of the astrophysical observations available up to date has 
been proposed on a purely phenomenological basis. It is, 
however, interesting to observe that the same model has a 
straightforward interpretation in the standard framework 
of a universe made out of two fluids, namely the matter 
term and the (dominant) dark energy. To see this, let us 
consider again Eq.Q. We may rewrite it as p = p\ + p 2 
with : 



For the model with (a, (3) = (3, 4), p\ scales with R as a 
matter -like term 5 , while pi as a quintessence - like fluid. 
Moreover, the matter -like term drives the evolution of 
the universe until R < b (i.e., z > zi,), after which the 
quintessence - like term starts dominating. Motivated by 
this analogy, we will refer to p± as p m and to p 2 as pq. Let 
us now define the two dimensionless density parameters : 



Q 



Pm(R=l) 
Pcrit 



pq(R = ±) 

Pcrit 



b a 



1 + b a 



(34) 



(35) 



Q m and £Iq may also be read as the density parameters 
of matter and quintessence, respectively. For a — 3, this 
gives : 

n m /n Q = b 3 ~ 0.37 -» b ~ 0.72 -v z b ~ 0.39 . 

On the other hand, using qo — —0.64 (as resulting from 
the angular size - redshift test using the Jackson data) 
and Eq. l(T5|) , we get : 



0A7^n m /n Q = (i + z b y 



0.32 



in qualitative good agreement with the previous estimate. 

Let us now evaluate the barotropic factors of these two 
fluids. To this aim, we may use Eq.@ to obtain: 



w m (R) 



w Q (R) 



(a-Z)R+(/3-Z)s 
3{R+s) 



3(R + s) + (a-P)s 



3(R + s) 



(36) 



(37) 



Note that, to use Eq.@, we have implicitly assumed that 
the two fluids do not interact so that the total pressure 
may be written as the sum of two single contributions. 
Actually it is quite easy to verify that : 



Pl {R) = Ab a R- p (R + s) f3 ~ a , 



(32) 



Ptot = Pm + Pq — w m p m + WqPq 

with ptot = wp and w given by Eq. (|10(l . 



P2 



(R) = A R a ~P {R + sf- a . 



(33) 



It is quite easy to see that : 



5 Here, as "matter - like" we mean here both radiation and matter. 
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It is worth noting that, for (a, f3) — (3,4), Eas. (|36() 
and l|37[) reduces to : 



w m (R) = 



3(R + . 



w Q (R) 



3(R + s) 



w m (R = 1) ~ , 



w Q (R=l) 



having used Eq.(J3J) and z s = 3454. The present day 
values of the barotropic factors for the two fluids only 
depend on the z s and reduce to that of matter and cos- 
mological constant, respectively, if z s ~ 10 -3 . 

Let us now concentrate on the quintessence - like term. 
In the standard framework, the quintessence fluid is gen- 
erated by a scalar field </> rolling down its potential V(<j>). 
These quantities are then related to the energy density 
Pq and the barotropic factor wq as follows || : 



PQ = \fr 



WQ = 



2V/j> 2 



(38) 



1 + 2V/<j) 2 



which, using Eq.lj^jl. can be solved 6 with respect to <j> and 
V to give : 



V Vn 



0.6 . 



T f fr 



FIG. 10: The scalar field potential for the model with (a, /3) = 
(3,4) and (z s ,q ) = (3454,-0.64). Here, (f> and Vo are the 
present day values of the scalar field and of the potential. 



the same as that of the standard cosmological constant. 
This result could also be expected from Eq. (|37|) . For the 
model with (a, (3) — (3, 4), the departures of wq from the 
the cosmological constant value w\ = — 1 are driven by 
the term (1 + R/s)^ 1 which is only slowly varying since 
R >> s over almost the full redshift range. It is worth 
stressing, however, that wq is never exactly equal to the 
cosmological constant value. This can also be seen con- 
sidering its derivative with respect to the redshift being : 



R" 



R h 



1 - w Q 

1 + 



dR' 

IF 



v 



WQ 



1 



wq 



WQ 



l+WQ 



PQ 



(39) 

Solving numerically this system allows to get the scalar 
field potential V(cf>) that we plot in Fig.^|for the model 
with (a,/3, z s ) — (3,4,3454). In this plot, we have fixed 
g = —0.64 as suggested by the best fit to the angular 
size-redshif test with the Jackson data. However, we 
have found that the shape of V/Vo vs </>/</>o does not de- 
pend on qo , while the deceleration parameter determines 
the present day values of both the scalar field (0o) and 
the potential (Vb). 

The most striking result is that, although a monoton- 
ically decreasing function of the scalar field as expected, 
the variation of the potential with respect to Vq is quite 
small, less than 2% over the full redshift range. In other 
words, the reconstructed scalar field potential is almost 



6 The reconstruction of the scalar field potential from the data is 
always possible provided that the right set of equations are used 
and the data are of sufficient quality to prevent strong degenera- 
cies in the potential reconstruction. An inter estin g example may 
be inferred from the equations presented in |27j where a model 
with varying scales and couplings is discussed. 



dw Q ((3-a)(l + z s ) 



dz 



{2 + z + z s f 



(40) 



For (a,f3) — (3,4) and whatever is the value of z s , this 
quantity does not identically vanish so that, even if the 
scalar field potential is slowly varying, the corresponding 
energy density may not be considered that of a cosmolog- 
ical constant (for which V((f>) does not depend at all on 
<f>). Therefore, some of the problems connected with the 
A phenomenology are not present for our parametriza- 
tion. On the other hand, a combined analysis of most 
of the available data gives the standard ACDM model as 
best fit Q so that it is not surprising that using a similar 
dataset individuates as best fit model among our class 
the one that best matches the ACDM one. 



VIII. CONCLUSIONS 

The concordance cosmological model assumes that 
pressureless cold dark matter and the dominant 
quintessence field drive the evolution of the universe lead- 
ing to the observed accelerated expansion. As an alter- 
native to this picture, unified dark energy models have 
been proposed where a single cosmic fluid acts both as 
dark matter and dark energy depending on the value of 
the energy density. 

Here we have presented a general class of models and 
tested it against the astrophysical observations available 
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up to date. The starting point is the assumption that 
the fluid energy density p depends on the scale factor R 
as shown in Eq.JIJ. The model is characterized by five 
parameters, but we have limited our attention to those 
models where the energy density smoothly interpolates 
among the three main phases of the universe evolution, 
i.e. a radiation dominated era followed by matter dom- 
ination and an asymptotic de Sitter state. This ansatz 
is motivated by the phenomenology we observe since ev- 
ery consistent picture of the universe evolution predicts 
the onset of the three different phases we have quoted 
above. The model has been tested against some astro- 
physical observations, namely the age of the universe, 
the SNela Hubble diagram and the angular size - rcdshift 
relation for compact radio structures. The succesful re- 
sults of these tests is a strong evidence of the reliability 
of this phenomenological approach. It is worth noting 
that our results may be interpreted both in the frame of 
UDE models and in the standard picture of a two fluid 
universe made of matter and dark energy. In this latter 
case, if a scalar field <j) is assumed to be the origin of 
this component, the interaction potential V(cj)) may be 
directly reconstructed from the observations without any 
a priori hypothesis on its form. 

Although the model has been shown to successfully fit 
the available observations, further analysis is still needed. 

First, there are other tests that can be used to bet- 
ter constrain the model parameters. In particular, we 
have not discussed the growth of perturbations which 
determines the CMBR anisotropy spectrum. Since the 
energy density p scales as the matter term in ACDM 
model during the structure formation era, perturbations 
should grow in the same manner so that the fit to the ob- 
served CMBR spectrum is likely to be successful. How- 
ever, this naive expectation has to be carefully checked 
since it depends on how we interpret the model. Ac- 
tually, this qualitative picture should indeed be valid if 
Eq.Q is meant as a phenomenological description of a 
two fluids scenario since, in this case, the quintessence - 
like fluid almost vanishes during the epoch of structure 
formation and therefore everyting works as in the usual 
scenario. On the other hand, if the UDE interpretation 
is preferred, one has to explicitely solved the perturba- 
tion equations by using explicitely Eq.Q for the matter 
enery density. We are, however, confident that the main 



results are not changed for the model with (a, (3) = (3, 4) 
since, in this case, Eq.QJ reduces to the usual expres- 
sion (p ax i?~ 3 ) in the structure formation epoch. Fi- 
nally, a different approach has to be considered if the 
interpretation of the model in the framework of modi- 
fied Friedmann equations. In this case, one should first 
obtain the corresponding modified Newtonian potential 
and then work out the consequences on the equations 
describing the growth of perturbation to finally obtain a 
coherent description of the structure formation process 
and of its imprint on the CMBR anisotropy spectrum. 

Secondly, in this analysis, we have held fixed the two 
slopes parameters (a, /?) to the values they must have to 
perfectly mimic the scaling of the energy density with 
the scale factor R during the radiation and matter dom- 
inated era. However, it is still possible that other values 
of (a,/3) fit the astrophysical data we have considered. 
It is thus interesting to repeat the same analysis per- 
formed here also varying these parameters. To this aim, 
however, more data are welcome since adding more pa- 
rameters may introduce strong degeneracies among some 
of them. In order to solve these problems, one should 
also consider the observations of large scale distribution 
of matter as measured by the two major ongoing galaxy 
surveys (the 2dFGRS 121 a nd the SDSS [H) or the data 
on the Lyman a forest 30]. 

As a final remark, we want to stress again the un- 
derlying philosophy of this work. Contrary to the usual 
approach of proposing a theory and then test it against 
observations, we have preferred to start from the phe- 
nomenology we observe to investigate what are the main 
features a theory should have to give a consistent and 
realistic picture of the universe. 
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